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[0001] This application claims priority to and the benefit of copending U.S. Provisional 
Patent Application Serial No. 60/21 1,461, filed on June 13, 2000, and entitled "A System and 
Method for Whole Genome Shotgun Assembly," the entire disclosure of which is hereby 
incorporated by reference herein* 
Government Support 

[0002] Work described herein was supported by Federal Grant Nos. 5 U54 HG021 52 and 
5 U54 HG02045, awarded by the National Institute for Health. The Government has certain 
rights in the invention. 
Technical Field 

[0003] The present invention relates generally to genetic information, and more 
specifically, to methods for the assembly of genomic sequences. 
Background of the Invention 

[0004] A problem encountered with assembly of genetic information from DNA 
sequences is that errors, e.g., base substitutions and indels, are introduced in the sequences 
during sequencing so that overlap between sequences can be difficult to detect. Another problem 
is that current sequencing programs can only resolve relatively short genome segments, i.e., the 
BAC clone that is roughly 200,000 nucleotides long. Another problem is the difficulty in 
assembling repeat-rich genomes from sequences, e.g., 25% of the human genome are repeated 
about 2 or 3 times within the genome, because the ambiguities are difficult to resolve. Because 
of the repeat problem, it is thought that sequencing repeat-rich genomes is not feasible. 
Summary of the Invention 

[0005] A novel approach to the assembly of genetic information has now been 
discovered, thus providing a method of assembly of genomic sequences. By practicing the 
disclosed invention, the skilled practitioner can now, e.g. , compare subsequences of the reads so 
that overlap can be detected despite sequence errors, significantly reduce the number of sequence 
comparisons necessary to detect overlaps, and resolve ambiguities due to the presence of 
repetitive regions in the genome by use of linking information. The advantages of the present 
invention include a feasible method of assembly of sequences that is efficient and trustworthy. 
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[0006] In one aspect, the invention features a method for assembly of a plurality of reads 
from a genomic region. The method includes the steps of providing a plurality of reads from a 
genomic region, indexing a plurality of read subsequences for each of the plurality of reads, each 
subsequence having an associated read with which it corresponds, extracting from the indexed 
subsequences a plurality of read pairs that have a predetermined number of subsequences in 
common, and merging the read pairs along a continuum. 

[0007] The method can include the step of providing a plurality of reads generated from 
sequencing both ends of a plurality of DNA segments, each read having associated linking 
information including an associated orientation relative to a read from an opposite end of the 
DNA segment, and an associated distance from the read on the opposite end of the DNA 
segment. The method can include the step of providing a plurality reads that are reverse 
complements of a plurality of reads provided by sequencing both ends of a plurality of DNA 
segments. The method can include the step of sorting the indexed read subsequences. The 
method can include the step of discarding read subsequences having more than a cutoff number 
of occurrences from the plurality of indexed subsequences. 

[0008] The method can include the step of indexing a plurality of read subsequences of a 
predetermined length for each of the plurality of reads. The predetermined length for each of the 
plurality of reads can be between about 12 and about 32 bases long. 

[0009] The method can include the step of indexing a plurality of subsequences for each 
read, the index including for each subsequence an associated read and an associated position on 
the read with which it corresponds. The method can include the step of performing alignments 
on the plurality of read pairs having a predetermined number of subsequences in common. The 
alignments can be performed on the plurality of read pairs having a predetermined number of 
subsequences in common and the associated position on the reads with which the subsequences 
correspond can be used to verify overlap. 

[0010] The method can include the step of using the linking information associated with 
the reads to confirm that the merged pairs are merged correctly. The method can include the step 
of using the associated linking information to an ambiguity in the merged reads. The method can 
include the step of identifying a repeat region and a set of unique regions. The method can 
include the step of linking pairs of unique regions using the linking information associated with 
the reads in the unique regions. The method can include the step of inserting the repeat region 
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between each linked pair of unique regions with which the repeat region corresponds. The 
method can include the step of merging linked pairs of unique regions using the linking 
information associated with the reads in the unique regions. 

[0011] In another aspect, the present invention features a method for assembly of merged 
reads from a genomic region. The method includes the steps of providing one or more sets of 
merged reads from a genomic region including a set of reads having associated linking 
information, using the associated linking information to identify an ambiguity in the merged 
reads, identifying a repeat region and a set of unique regions, and linking pairs of unique regions 
using the linking information associated with the reads in the unique regions. 

[0012] The method can include the step of inserting the repeat region between each 
linked pair of unique regions with which the repeat region corresponds. The method can include 
the step of 

[0013] merging linked pairs of unique regions using the linking information associated 
with the reads in the unique regions. 

[0014] In yet another aspect, the present invention features an article of manufacture 
having computer-readable program means embodied thereon for assembly of a plurality of reads 
from a genomic region. The article includes computer-readable program means for providing a 
plurality of reads from a genomic region, computer-readable program means for indexing a 
plurality of read subsequences for each of the plurality of reads, each subsequence having an 
associated read with which it corresponds, computer-readable program means for extracting from 
the indexed subsequences a plurality of read pairs that have a predetermined number of 
subsequences in common, and computer-readable program means for merging the read pairs 
along a continuum. 

[0015] In yet another aspect, the present invention features an article of manufacture 
having computer-readable program means embodied thereon for assembly of merged reads from 
a genomic region. The article includes computer-readable program means for providing one or 
more sets of merged reads from a genomic region including a set of reads having associated 
linking information, computer-readable program means for using the associated linking 
information to identify one or more ambiguities in the merged reads, computer-readable program 
means for identifying a repeat region and a set of unique regions, and computer-readable 
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program means for linking pairs of unique regions using the linking information associated with 
the reads in the unique regions. 

[0016] In short, the invention provides the art with a heretofore unappreciated method for 
assembling genomic sequences that allows for sequencing errors, reduces the number of read 
comparisons, and resolves ambiguities due to repetitive regions in the genome. 

[0017] These and other objects of the present invention, along with advantages and 
features of the invention disclosed herein, will be made more apparent from the description, 
drawings and claims that follow. 
Brief Description of the Drawings 

[0018] The invention is pointed out with particularity in the appended claims. The 
drawings are not necessarily to scale and simplified for the sake of clarity and brevity, emphasis 
instead generally being placed upon illustrating the principles of the invention. The advantages 
of the invention can be better understood by reference to the description taken in conjunction 
with the accompanying drawings, in which: 

[0019] Figures 1A-D are conceptual diagrams illustrating an exemplary embodiment of a 
method of providing a plurality of reads from a genomic region in accordance with the present 
invention; 

[0020] Figure 2 is an illustration of a exemplary embodiment of a method of indexing a 
plurality of read subsequences for each of a plurality of reads, each subsequence having an 
associated read with which it corresponds, in accordance with the present invention; 

[0021] Figure 3 is an illustration of an exemplary embodiment of a method of extracting 
from indexed subsequences a plurality of reads pairs that have a predetermined number of 
subsequences in common, in accordance with the present invention; 

[0022] Figure 4 is a illustration of an exemplary embodiment of a method of merging 
read pairs along a continuum in accordance with the present invention; 

[0023] Figures is a illustration of an exemplary set of merged reads having a repeat 
region and a set of four unique regions in accordance with the present invention; 

[0024] Figure 6 is an illustration of an exemplary embodiment of a method of using 
linking information to identify an ambiguity in the merged reads in accordance with the present 
invention; 
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[0025] Figure 7 is an illustration an exemplary embodiment of a method of identifying a 
repeat region and a set of unique regions, in accordance with the present invention; 

[0026] Figure 8 is an illustration of an exemplary embodiment of a method of linking 
pairs of unique regions using the linking information associated with the reads in the unique 
regions, and inserting the repeat region between each linked pair of unique regions with which 
the repeat region corresponds, in accordance with the present invention; and 

[0027] Figure 9 is an illustration of an exemplary method of merging linked pairs of 
unique regions using the linking information associated with the reads in the unique regions in 
accordance with the present invention. 
Detailed Description 

[0028] In its broadest aspects, the present invention is directed to the reconstruction or 
assembly of an unknown source sequence of a set of nucleotides or bases {A,C,G,T} from a 
genomic region, given a plurality of random sequences or "reads" from the sequence. Reads 
typically are about 500 nucleotides long, but can be in the range from about 200 to about 1000 
nucleotides long. The methods of the present invention contemplate the indexing of read 
subsequences for each read. Read pairs are then extracted from the index that have a number of 
subsequences in common. These read pairs are then merged along a continuum. This method 
allows for or tolerates sequencing errors because subsequences of each read are compared, as 
opposed to the entire read sequence. This method obviates the need to compare every read to 
every other read, reducing an N 2 number of calculations to a linear number as the number of 
reads increases. If the DNA is several billion nucleotides long, transformation of a polynomial 
N problem to a linear problem using the methods of the present invention represents a 
significant improvement in efficiency and speed. 

[0029] The methods of the present invention also contemplate the assembly of read pairs 
obtained from both ends of DNA segments. Each read in the pair has associated with it 
information including the approximate distance between the reads and the relative orientation of 
the reads. This information can be used to identify repetitive regions in the DNA and resolve 
ambiguities in the merged reads. 

[0030] One aspect of the present invention is a method for assembly of a plurality of 
reads from a genomic region. The method generally includes the steps of: providing a plurality 
of reads from a genomic region; indexing a plurality of read subsequences for each of the 
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plurality of reads, each subsequence having an associated read with which it corresponds; 
extracting from the indexed subsequences a plurality of read pairs that have a predetermined 
number of subsequences in common; and merging the read pairs along a continuum. 

[0031] The genomic region can be from any source of DNA including human DNA. The 
genomic region can be the entire genome or a portion of the genome. Generally, the greater the 
depth of read overlap, the greater the confidence in the assembly. If the genomic region is length 
L, preferably enough DNA copies are provided so that the sum of the lengths of the reads is 
between about 2L and 20L. More preferably, enough DNA copies are provided so that the sum 
of the lengths of the reads is about 10L. The plurality of reads from a genomic region can be 
provided by known methods, e.g., by shotgun sequencing using known methods. 

[0032] An embodiment of this method includes the step of providing a plurality of reads 
generated from sequencing both ends of a plurality of DNA segments, each read having 
associated linking information. This linking information includes an associated orientation 
relative to a read from an opposite end of the DNA segment, and an associated distance from the 
read on the opposite end of the DNA segment. 

[0033] Figures 1 A-D are conceptual diagrams illustrating an exemplary embodiment of a 
method of providing a plurality of reads from DNA in accordance with the present invention. 
Figure 1A depicts 5 copies of a DNA sequence 1. The DNA can be any DNA, e.g., human 
DNA, which is roughly 3 million base pairs long. The DNA is broken up into random segments 
4 as shown in Figure IB using methods known in the art. The segments 4 then are inserted into 
appropriate vectors, e.g., plasmids or cosmids, 8 using known methods as shown in Figure 1C. 
Then, referring to Figure ID, using known methods, e.g., gel electrophoresis, typically a 
sequence of approximately 500 nucleotides is read from each segment end 12, 16. Thus, for 
each segment 4, a pair of reads 20, 24, is obtained with a known orientation relative to each other 
and approximate distance d from each other. That is, it is known that: read 20 has the sequence 
ACGTA . . . CCGTTAT; the "T" end of the read 20 is oriented toward read 24; the read 24 has 
the sequence AGT . . . ATACGGAT; the "T" end of read 24 is oriented towards read 20; and the 
approximate distance, d, between read 20 and 24 because we know the length of segment 4 
inserted into vector 8. The reads are sequenced from the ends of segment 4, one from the 
forward strand of DNA and one from the cDNA strand. 
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[0034] The segments can be of varying length. For example, the segments can all be 
relatively short, e.g., 2,000-3,000 nucleotides long, or relatively long, e.g., about 150,000 
nucleotides long. Alternatively, the segments can include groups of varying segment lengths, 
e.g., in an embodiment of the present invention, approximately one-tenth of the segments used to 
generate the reads are about 40,000 nucleotides long and approximately nine-tenths of the 
segments are about 4,000 nucleotides long. 

[0035] The method preferably includes the step of providing a plurality of reads that are 
reverse complements of the plurality of reads provided by sequencing both ends of a plurality of 
DNA segments. These reverse complement reads can be provided by reversing the nucleotide 
sequence order and substituting each base with its Watson-Crick complement: AoT and C <-» 
G. Each of these reverse complement reads preferably are associated with the linking 
information discussed above including the orientation of the read relative to the read on the other 
end of the segment, and the approximate distance from the read to the read on the other end of 
the segment. 

[0036] Providing the reverse complement reads is preferred because the relative direction 
of each read, or pair of reads obtained from reading off both ends of a DNA segment, is not 
known relative to the other reads. Without these reverse complement reads, only about half of 
the read overlaps would be detected when merging the read subsequences. For example, 
consider a read that includes the subsequence CCGAATGA. Reads including the sequence 
CCGAATGA can readily be detected as possibly overlapping. However, reads including the 
reverse complement sequence TCATTCGG would not be detected. The plurality of reads 
provided can also include reads without any associated link information and/or reads taken from 
known portions of the genomic region. 

[0037] The method includes the step of indexing a plurality of read subsequences for 
each of the plurality of reads, each subsequence in the index is associated with its corresponding 
read. Each subsequence in the index also can have a "forward" or "reverse" direction associated 
with it. 

[0038] Figure 2 is an illustration of an exemplary method of indexing read subsequences 
in accordance with the present invention. Figure 2 depicts a Table of Reads and a Table of 
Subsequences. The Table of Reads includes each read sequence and a read number assigned 
each read for identification purposes. For the purposes of this illustration, only 6 reads are 
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shown and the reads are only 8 nucleotides long, however, in practice typically there are many 
more reads, each about 500 nucleotides long. This table also includes an indication of whether 
each sequence is from the forward or reverse complement read. From the Table of Reads, an 
index of subsequences is created, as illustrated in the Table of Subsequences. The Table of 
Subsequences includes every subsequence of each read. 

[0039] Each subsequence is associated with a read number and a read direction 
corresponding to whether the subsequence corresponds to a read provided from reads generated 
from sequencing both ends of a DNA segment (F), or a reverse complement of these reads (R). 
For example, read number 1 has the sequence CCTGCGCC. It has 5 subsequences that are 4 
bases long. Its 4-mers are: CCTG, CTGC, TGCG, GCGC and CGCC. These subsequences are 
depicted as the first 4 entries in the Table of Subsequences associated with its read number and 
read direction. The reverse complement of read 1 is GGCGCAGG. It's 4-mers are GGCG, 
GCGC, CGCA, GCAG and CAGG. These subsequences are not depicted in Figure 2 for the 
sake of brevity, but would also be included in the index in the practice of this embodiment. 

[0040] For the purposes of illustration, the subsequences in Figure 2 are 4 nucleotides 
long. In practice, the subsequences can be of any predetermined length (k-mers) or of random 
lengths. If a predetermined length is used, preferably its between about 12 and about 32 
nucleotides long, more preferably between about 20 and about 27 nucleotides long, and most 
preferably about 24 nucleotides long. These lengths allow or tolerate for sequencing errors that 
occur in about 2% of the nucleotides. For example, consider a first read of 500 nucleotides that 
has errors at nucleotide numbers 2, 50, 99, 157, 206, 240, 359, 361, 399 and 41 1, and an 
overlapping read of 520 nucleotides with errors at nucleotides corresponding to the first read's 
nucleotide numbers 8, 15, 120, 165, 199, 280, 310, 330, 450 and 458. There are several regions 
where 24-mers can be error-free and thus overlap between 24-mers readily can be detected, e.g., 
between reads 16-49, 51-98, 1 19-156 and 166-198. 

[0041] The Table of Subsequences depicted in Figure 2 can include additional 
information. For example, the Table of Subsequences can include information about the position 
of the subsequence along the read. As discussed below, this information can be useful when 
performing alignments on read pairs so that overlap can be verified. 

[0042] The indexed subsequences can be sorted as shown in Figure 3. The sort can be 
by nucleotide. The sort can be used to create an index where identical subsequences are adjacent 
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to one another in the index so that read pairs having subsequences in common readily can be 
identified. For example, Figure 3 shows us that subsequence AGCC is common to reads 3 and 
61, and subsequence CCCT is to reads 3, 50 and 61. 

[0043] The method can further include the step of discarding read subsequences having 
more than a cutoff number of occurrences, from the indexed subsequences. This step can be 
used, e.g., to remove from the index subsequences that are repeated so many times that it is not 
efficient to include them for the purposes of detecting read overlaps. For example, in the human 
genome, the nucleotide sequence for the amino acid Alu is repeated as many as a million times 
throughout the genome, and so it is not efficient to consider the nucleotide codign for Alu 
subsequences when attempting to detect read overlaps. The cutoff number can be determined by 
the initial redundancy of sequencing that was used. For instance, if the total length of the reads 
is 10L, where L is the length of the genome in nucleotides, the initial redundancy of sequencing 
is 10, and therefore an average region of the genome is expected to be represented 10 times by 
reads. Therefore, a subsequence that occurs considerably more than 10 times is likely to be 
repetitive. A subsequence that occurs 100 times is likely to be so repetitive that it is not efficient 
to include it for the purposes of detecting overlaps, and such subsequences can be discarded. 
The best cutoff number is the largest cutoff number that allows for an acceptably low number of 
pairwise comparisons. 

[0044] A Table of Sorted Sequences, such as that depicted in Figure 3, can also include a 
"Remove Flag" column (not shown) that identifies subsequences that occur more than a cutoff 
number of times. This step can be used to remove highly repetitive subsequences from 
consideration to increase the efficiency of the merging step, since these subsequences likely will 
falsely indicate overlap between the reads in which they appear. 

[0045] The method further includes the step of extracting from the indexed subsequences 
a plurality of read pairs that have a predetermined number of subsequences in common. The 
common subsequence or subsequences indicate that the read pair has at least some region or 
regions of overlap. This method reduces the number of read comparisons that have to be 
performed to merge the pairs from the number of reads squared (N 2 ) to the number of pairs that 
have one or more common subsequence. 

[0046] The predetermined number of common subsequences can be 1 or higher. 
Preferably, the predetermined number is 1 . However, if a collection of reads having a large 
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percentage of errors were to be assembled, the subsequence length can be small and the 
predetermined number greater than 1 so that overlaps could be detected despite the errors. 
Generally, the higher the predetermined number, the greater the chance that a region of overlap 
can be missed, and the smaller the number of extracted read pairs. 

[0047] Figure 3 is an illustration of an exemplary method of extracting read pairs having 
common subsequences by building a Table of Read Pairs with Common Subsequences from a 
sorted Table of Sorted Subsequences in accordance with the present invention. The Table of 
Sorted Subsequences is a sorted index of the Table of Subsequences depicted in Figure 2. The 
Table of Read Pairs with Common Subsequences is generated from the Table of Sorted 
Sequences by extracting a list of read pairs having at least 1 subsequence in common. Each read 
pair is associated with the number of common subsequences indicated in parenthesis. 

[0048] Figure 3 illustrates that this method of indexing the subsequences and extracting 
the read pairs having common subsequences significantly reduces the number of read 
comparisons for the merge from the number of reads squared (N 2 ), to the number of read pairs 
that have one or more common subsequences. For example, using the method of the present 
invention, read 1 will be compared with reads 3 and 50, but not with 61 or 14, whereas with an 
N comparison method, read 1 would be compared with every other read, including 14 or 61 . 
Similarly, because read 105 has no subsequences in common with reads 1, 3, 14, 50 and 61, 
these comparisons are not made. This method can include sorting the table of subsequences by 
nucleotide. 

[0049] The method further includes the step of merging the read pairs along a continuum. 
The pairs extracted from the index are compared; if they overlap, they are aligned and merged 
along a continuum. This step can include the step of performing alignments on the extracted 
read pairs to confirm that the reads overlap. Methods of performing alignments are known. 

[0050] Figure 4 is an illustration of an exemplary method of merging read pairs along a 
continuum. This merge can be performed by choosing one read pair, e.g., read pair (3,61) from 
Figure 3. A direction is arbitrarily chosen for the first alignment, in this case, it is decided to 
represent the reads from left to right. Pair (3,61) is aligned along a continuum running from left 
to right. Next, a read pair can be chosen that includes one of the reads in the pair already 
aligned. Pair (3,50) can be chosen, e.g., and read 50 aligned with read 3 in Figure 4. The merge 
can be built pair by pair as shown in Figure 4. It should be noted that read 14, shown in Figure 
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4, would not be aligned from the data shown in Figure 3 because it does not share a common 
subsequence with the other reads shown. However, as the reads are merged with data not shown 
in Figure 4, other pairs can be aligned until a read is added that has a subsequence in common 
with read 14, and read 14 would be merged. From the merged reads, the resolved sequence 
ATAGCCCTGCGCCTATCG, indicated below the continuum line in Figure 4, can be obtained. 

[0051] Alignments optionally can use the associated position of the subsequences on the 
reads to confirm overlap. For example, consider the two sequences: GATCCCATGCGCA and 
ATAGCCCTATGAT. These sequences share common subsequences GAT and CCC. We know 
from the associated position information that CCC begins at base 4 for the former and base 5 for 
the latter, and the GAT subsequence begins at base 1 for the former and base 1 1 for the latter. 
This position information indicates that there is no overlap between these two sequences since 
the position of the common subsequences does not allow for alignment. Consider now the first 
sequence above and the sequence CCCATGCGC AT AT. We know that the common 
subsequence CCC begins at base 4 for the former and base 1 for the latter, and the common 
subsequence CGC begins at base 9 for the former and base 6 for the latter. This position 
information indicates that there is overlap between these two sequences. Comparing the rest of 
the intervening subsequences confirms the overlapping region CCCATGCGCA. 

[0052] The method can further include the step of using the linking information 
associated with the reads to confirm that the merged pairs are merged correctly. The linking 
information can include the relative orientation of the reads and the approximate distances 
between the reads. This information can be used to confirm the merged reads, by checking for 
some or all of the reads, whether they have been merged such that the relative orientation of the 
reads and the distance between the reads are consistent with the linking information. 

[0053] The merge can contain ambiguities as encircled in Figure 5, where it appears that 
there is a divergence in the sequence. These ambiguities are due to repeat regions in the DNA, 
and are discussed below. 

[0054] The merge can also contain discontinuities (not shown). That is, there can be 
gaps in the merge where regions of overlap were not identified and thus the merge can yield two 
or more merged regions of pairs. From the merged pairs alone, the relative orientation of the 
merged regions along the continuum is unknown. However, as discussed below, the relative 
orientation of these regions can be determined using linking information. 
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[0055] As noted above, Figure 5 depicts a read merge that contains encircled ambiguities 
where the sequence appears to diverge. Figure 5 is a illustration of an exemplary set of merged 
reads, represented by overlapping lines, having a repeat region R and a set of four unique regions 

A, B, C, and D. This merge is ambiguous because it is not known how the regions should are 
linked because of the presence of repeat region R. That is, it is not known whether regions A, C, 

B, D should be combined such that regions A-R-B and C-R-D are linked in the genome, or A-R- 
D and regions C-R-B are linked in the genome. Figure 5 depicts a merge having 4 unique 
regions A, B, C, D for the sake of illustration, however, in practice merges can have a large 
number of unique regions when repetitive region R is repeated a large number of times in the 
genome. How these regions are linked can be difficult, if not impossible, to determine with read 
subsequences alone. 

[0056] Another aspect of the present invention is a method for assembly of merged reads 
from a genomic region. This method addresses the difficulty of ordering ambiguous merges. 
The method generally includes the following steps: providing one or more sets of merged reads 
from a genomic region, the merged reads including a plurality of reads generated from 
sequencing both ends of a plurality of DNA segments, each read having associated linking 
information including an associated orientation relative to a read from an opposite end of the 
DNA segment, and an associated distance from the read on the opposite end of the DNA 
segment; identifying one or more ambiguities in the merged reads where the merged reads 
diverge; using the associated linking information to identify a repeat region and a plurality of 
unique regions; and linking pairs of unique regions using the linking information associated with 
the reads in the unique regions. 

[0057] The one or more sets of merged reads provided can include merges provided 
using the methods of the present invention as described above. Additionally or alternatively, the 
set of merges can include merges built using other known methods such as those created by 
performing N alignments. At least some of the merged reads include a plurality of reads 
generated from sequencing both ends of a DNA segment as described above. These reads each 
have associated linking information including an associated orientation relative to a read from an 
opposite end of the DNA segment, and an associated distance from the read on the opposite end 
of the DNA segment. Optionally, the merged reads can include reads that are not generated from 
sequencing both ends of a DNA segments. Optionally, the merged reads can include sequences 
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obtained using other sequencing techniques known in the art, e.g., nucleotide tagging and 
sequencing, provided that linking information in known for at least some of the "reads" 
embedded in the sequence. 

[0058] The method includes the step of using the associated linking information to 
identify one or more ambiguities in the merged reads. An ambiguity is detected when 
overlapping reads have read pairs that are contained in two distinct merged regions which, 
according to the linking information, should also overlap, but do not. 

[0059] Figure 6 is an illustration of an exemplary method of using the associated linking 
information to identify one or more ambiguities in the merged reads, in accordance with the 
present invention. Figure 6 depicts an exemplary set of merged reads ri, r 2 , r 3 i, r 4 3, r 50 and r 6 . 
The ambiguity can be detected by choosing two reads in the same set of merged reads, e.g., reads 
ri and rso. The reads in Figure 6 are indicated as r„, where n is the read number. In this 
example, the reads were generated from sequencing both ends of both strands of a DNA segment 
and so approximate distance di,5 between r\ and rs and distance dso,83 between rso and rs3, and 
relative orientation of the reads, as indicated by the arrow heads are known. Reads rs and rs3 are 
in unique merged regions B and D. This information tells us that B and D should overlap, 
however, a comparison of the sequences of B and D show that there is no overlap. This indicates 
that there is an ambiguity in the reads. A similar analysis of linking information between, e.g., 
r4o and r 3 i, and r 75 and r$, tells us that A and C should overlap, however a comparison of these 
sequences show that there is no overlap, indicating that there is an ambiguity in the reads. 

[0060] The method further comprises the step of identifying a repeat region and a set of 
unique regions. Once the ambiguity is identified, a review of the sequences of the merged reads 
about the ambiguous regions yields the identity of the repeat and unique regions. Generally, each 
"arm" or branch about the ambiguous regions is a unique region, and the region between the 
ambiguous regions is the repeat region. 

[0061] Figure 7 is an illustration of an exemplary method of identifying a repeat region 
R and a set of unique regions A, B, C, D, in accordance with the present invention. By distilling 
the sequences of the merged reads, as shown in Figure 7, the portion of the merged reads that is 
not branched can be identified as repeat region R, and the arms or branches are unique regions 
A, B, C, D. Only two unique regions are shown on each side of repeat region R for the sake of 
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illustration. However, in practice a repeat region can be flanked by more unique regions 
depending on the number of times the repeat region appears in the genome. 

[0062] The correct assembly of the regions is not known from the information depicted 
in Figure 7. That is, it is not known if regions A and B flank one copy of repeat region R and 
regions C and D the other, or whether unique regions A and D flank one copy of repetitive 
region R and C and B the other. 

[0063] The method further comprises the step of linking pairs of unique regions using the 
linking information associated with the reads in the unique regions. Using linking information 
between reads in the unique regions, the correct assembly of the unique regions can be 
determined. After the unique regions are linked, if the repeat region was removed from between 
the unique regions, the method can comprise the step of inserting the repeat region between each 
linked pair of unique regions with which the repeat region corresponds. 

[0064] Figure 8 is an illustration of an exemplary method of linking pairs of unique 
regions using the linking information associated with the reads in the unique regions, and 
inserting the repeat region between each linked pair of unique regions with which the repeat 
region corresponds, in accordance with the present invention. The correct assembly of the 
unique regions can determined by searching or otherwise identifying reads in the unique regions 
with linking information to reads also in the linking regions. For example, Figure 8 depicts a 
link between reads 111 and 62 (rm, r^), with a known orientation relative to each other and a 
known distance between reads, dm^. It also shows reads 320 and 305 (r 32 o, r 30 s) 5 with a known 
orientation relative to each other and a known distance between the reads d 3 203os. From this 
linking information, it can be determined that regions A and D flank one copy of repeat region 
R, and regions B and C flank a second copy of repeat region R. Once the linked pairs are 
identified, the method can optionally include inserting the repeat region between the linked pairs 
it corresponds to as shown in Figure 8. What is not known from the information depicted in 
Figure 8 is how the linked pairs of unique regions are oriented relative to each other and at what 
distance. 

[0065] The method of the present invention can include the step of merging linked pairs 
of unique regions using the linking information associated with the reads in the unique regions. 
The linked pairs of unique regions to be merged can be from several sets of merged reads having 
ambiguities. 
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[0066] Figure 9 is an illustration of an exemplary method of merging linked pairs of 
unique regions using the linking information associated with the reads in the unique regions in 
accordance with the present invention. In Figure 9, unique regions D and G are linked by reads 
in these regions that are generated from both ends of a DNA segment having an approximate 
distance and known orientation relative to each other. This information indicates the relative 
distance and orientation of unique regions A-D and unique regions G-H. Merge 1 is possible 
because G-H can be placed on a continuum with A-D without causing a divergence in the 
sequence. The dashed lines indicate a repetitive region. 

[0067] In Merge 2, unique regions D and O are linked by reads in these regions as well 
This information indicates the relative distance and orientation of A-D-G-H and O-P. Merge 2 
is possible because A-D-G-H can be merged without causing a divergence in the sequence. 

[0068] In Merge 3, unique regions A and R are linked by reads in these regions. 
However, this time the distance associated with the link is not consistent with the presence of A- 
D-O-G-H-P, as shown by the circled region, and Merge 3 is not performed. This inconsistency 
can occur, e.g., if region A or R is a repetitive region, leading to inconsistent linking information. 
At this point the merging stops so that the ordering of non-repetitive regions is not compromised. 

[0069] Practice of the invention will be still more fully understood from the following 
example, which are presented herein for illustration only and should not be construed as limiting 
the invention in any way. 

Example: Simulated assembly of Chromosome 22 of the Human Genome 

[0070] The method was used to assemble shotgun reads from the chromosome 22 of the 
human genome. Simulated reads were created from the real sequence of chromosome 22. All of 
the reads came from earmuff links. Approximately 350,000 random segments were taken from 
approximately 10 copies of chromosome 22. Approximately one-third of these segments were 
about 10,000 bases long and about two-thirds of the segments were about 3,000 bases long. The 
former were simulating the insertion into long plasmid vectors and the latter into short plasmid 
vectors, and the approximately the first and last 500 bases from each end of the segments were 
sequenced. Sequencing errors were distributed uniformly at random and were present in an 
average error rate of 0.5% for mutations, 0.05% or insertions and 0.05% for deletions. The 
distance between the reads on the segments and the relative orientation of the reads were 
recorded. The earmuff link distances had an average error of 5%. Reverse complement reads 
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were provided that were reverse complements of the reads provided by sequencing both ends of 
the DNA segments. 

[0071] An index of 700,000,000 subsequences 24 bases long was generated for each of 
the reads, each having an associated read number. The index was sorted according by 
subsequence, and subsequences that had more than 100 occurrences were discarded. From this 
index, approximately 28,000,000 pairs (approximately 40 times the number of reads) were 
extracted that had at least 1 subsequence in common. 

[0072] These read pairs were merged along a continuum to provide 4,200 sets of merged 
reads. Ambiguities were identified in 1,300 sets of merged reads using linking information. The 
repeat regions and unique regions were identified. The unique regions were linked using linking 
information associated with reads in the unique regions. The linked pairs were merged using 
linking information associated with the reads in the unique regions. The repetitive regions were 
then inserted between each linked pair of unique regions with which the repeat corresponds. 
This method was processed by a single Alpha processor with 2Gb memory and 50Gb storage and 
took about 15 hours to complete. 
Equivalents 

[0073] The present invention can be implemented as one or more computer-readable 
software programs embodied on or in one or more articles of manufacture. The article of 
manufacture can be, for example, any one or combination of a floppy disk, a hard disk, hard-disk 
drive, a CD-ROM, a DVD-ROM, a flash memory card, an EEPOM, an EPROM, a PROM, a 
RAM, a ROM, or a magnetic tape. In general, any standard or proprietary, programming or 
interpretive language can be used to produce the computer-readable software programs. 
Examples of such languages include C, C++, Pascal, JAVA, BASIC, Visual Basic, and Visual 
C++. The software programs may be stored on or in one or more articles of manufacture as 
source code, object code, interpretive code, or executable code. 

[0074] While the invention has been shown and described with reference to specific 
preferred embodiments, it should be understood by those skilled in the art that various changes in 
form and detail may be made therein without departing from the spirit and scope of the invention 
as defined by the following claims. 
What is claimed is: 
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